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Abstract. In this paper we study the effects of rotation in low-mass, low-metallicity RGB stars. We present the first evolutionary 
models taking into account self-consistently the latest prescriptions for the transport of angular momentum by meridional 
circulation and shear turbulence in stellar interiors as well as the associated mixing processes for chemicals computed from 
the ZAMS to the upper RGB. We discuss in details the uncertainties associated with the physical description of the rotational 
mixing and study carefully their efi^ects on the rotation profile, diffusion coefficients, structural evolution, lifetimes and chemical 
signatures at the stellar surface. We focus in particular on the various assumptions concerning the rotation law in the convective 
envelope, the initial rotation velocity distribution, the presence of /j-gradients and the treatment of the horizontal and vertical 
turbulence. 

This exploration leads to two main conclusions : (1) After the completion of the first dredge-up, the degree of differential 
rotation (and hence mixing) is maximised in the case of a differentially rotating convective envelope (i.e., jcE(r) = est), as 
anticipated in previous studies. (2) Even with this assumption, and contrary to some previous claims, the present treatment for 
the evolution of the rotation profile and associated meridional circulation and shear turbulence does not lead to enough mixing of 
chemicals to explain the abundance anomalies in low-metallicity field and globular cluster RGB stars observed around the bump 
luminosity. This study raises questions that need to be addressed in a near future. These include for example the interaction 
between rotation and convection and the trigger of additional hydrodynamical instabilities. 
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1. Abundance anomalies in RGB stars 

The standard theory of stellar evolution' predicts that the sur- 
face chemical composition of low-mass stars is modified on 
the way to the red giant branch (RGB) during the so-called 
first dredge-up (hereafter 1st DUP; Iben ll965l . There, the ex- 
panding stellar convective envelope (hereafter CE) deepens in 
mass, leading to the dilution of the surface material within re- 
gions that have undergone partial hydrogen burning on the ear- 
lier main sequence (hereafter MS). Qualitatively, this leads to 
the decrease of the surface abundances of the fragile LiBeB el- 
ements and of '^C, while those of ^He, '^C and '''N increase. 
Abundances of O and heavier elements remain essentially un- 
changed. Quantitatively, these abundance variations depend 
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' By this we refer to the modelling of non-rotating, non-magnetic 
stars, in which convection and atomic diffusion are the only transport 
processes considered. 



on the stellar mass and metallicity (e.g., Sweigart, Greggio 
& Renzini,[nH3 Charbonnel FD^ Boothroyd & Sackmann 
\T$99 ) . After the 1 st DUP, the CE withdraws while the hydrogen 
burning shell (hereafter HBS) moves outward in mass. Within 
the standard framework no more variations of the surface abun- 
dance pattern are expected until the star reaches the asymptotic 
giant branch. 

Observations sampling the evolution from the turn-off to 
the base of the RGB in open clusters and in the galactic field 
stars have validated these predicted surface abundances vari- 
ations up to the completion of the 1st DUP^ (e.g. Gratton et 
al. I2000> . However observational evidence have accumulated 
of a second and distinct mixing episode which is not predicted 



- One has of course to take into account possible variations of the 
surface abundance of lithium occurring in some cases already on the 
MS. This discussion is however out of the scope of this paper (see 
e.g. Charbonnel, Deliyannis & Pinsonneault I2000l and Palacios et al. 
|2003, paper 1\ . 
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by standard models and which occurs in low-mass stars after 
the end of the 1st DUP, and more precisely at the RGB bump. 

The determination of the carbon isotopic ratio '^C/'^C 
(hereafter CIR) for RGB stars in open clusters with various 
turn-off masses (Gilrov [T989> provided the first pertinent clue 
on this process. It was indeed shown that bright RGB stars with 
initial masses lower than ~ 2 - 2.5 M© exhibit CIR consider- 
ably lower than predicted by standard models after the 1 st DUR 
Thanks to data collected in stars sampling the RGB of M67 
(Gilroy & Brown n"991> . it clearly appeared that observations 
deviated from standard predictions just at the so-called RGB 
bump (Charbonnel l 1 994l l . The Hipparcos parallaxes allowed to 
precisely determine the evolutionary phase of large samples of 
field stars with known CIR. These stars were found to behave 
similarly as those in M67, e.g. presented unpredicted low CIR 
appearing at the RGB bump luminosity (Charbonnel, Brown & 
Wallerstein [T998l Gratton et al.L2000l. 

On the other hand the region around the bump has also been 
probed for two globular clusters (GCs). In NGC 6528 and M4 
again, the CIR drops below the 1st DUP standard predictions 
just at the RGB bump JShetrone 2()03all2003bl l. Moreover, all 
the brightest RGB stars observed so far in globular clusters ex- 
hibit CIR close to the equilibrium value of the CN cycle. 

During this second mixing episode surface abundances of 
other chemical elements are also affected both in field and GCs 
giants : Li decreases at the RGB bump (Pilachowski, Sneden 
& Booth [T9931 Grundahl et al. l2002l . C decreases while N in- 
creases for RGB stars brighter than the bump (Gratton et al. 
120001 Bellman et al. 1200 II and references therein), confirming 
the envelope pollution by CN processing. In the case of GCs, 
the picture is however blurred by the probable non-negligible 
dispersion of the initial [C/Fe]. As far as lithium, carbon iso- 
topes and nitrogen are concerned, the abundance variations on 
the upper RGB have similar amplitudes in field and globular 
cluster giants (Smith & Martell 12003 ). The finding that the so- 
called super Li-rich giants (Wallerstein & Sneden 1982) all lie 
either at the RGB bump or on the early-AGB (Charbonnel & 
Balachandran I2000l . certainly indicates the occurrence of an 
extra-mixing episode at these evolutionary points. The trigger 
of this mixing episode has been suggested to be of external na- 
ture (Denissenkov & Herwig, 2004| l, but is much likely related 
to the aforementioned second mixing episode, which would 
start with a Li-rich phase as proposed by Palacios et al. (1200 U . 

For more than a decade it has been known that in addition 
to the elements discussed previously, O, Na, Mg and Al also 
show variations in GC red giants (Kraft et al. 119931 Ivans et 
al. 119991 Ramirez & Cohen EUU^ . As in the case of lighter 
nuclei, an in situ mixing mechanism was frequently invoked to 
explain these abundance anomalies, and in particular the O-Na 
anti-correlation. For a long time this specific pattern could 
only be observed in the brightest GC RGB stars. However, 
O and Na abundances have recently been determined with 
8-lOm class telescopes in lower RGB and in turn-off stars 
for a couple of GCs (for recent reviews see Sneden 2005 and 
Charbonnel 2005 and references therein) revealing exactly 
the same O-Na anti-correlation as in bright giants. This 
result is crucial. Indeed, the NeNa-cycle does not operate in 
MS low-mass stars, as the involved reactions require high 



temperatures that can only be reached on the RGB. The 
existence of the same O-Na anti-correlation on the MS and 
on the RGB in these clusters thus proves that this pattern 
does not result from self-enrichment. The recent determination 
of oxygen isotopic ratios in a few RGB stars with low CIR 
(Balachandran & Carr I2003> reinforces this result. These 
objects indeed present high '^O/'^O and '*0/"*0 ratios, in 
agreement with extensive CN-processing but no dredge-up 
of ON-cycle material. The O-Na anti-correlation is gener- 
ally assumed to be of primordial origin, even though it has 
been proved difficult to find stellar candidates able to produce it 
dDecressin & Charbonnel, 2005|Denissenkov & Herwig, 2003| l. 
Let us finally mention the peculiar case of Ml 3 bright giants, 
where O, Na, Mg and Al abundances appear to vary with 
luminosity. In this cluster the observed O-Na anti-correlation 
could thus be the result of a superimposition of self-enrichment 
with a primordial pattern (?). 

In short, observations provide definitive clues on an addi- 
tional mixing episode occurring in low-mass stars after the end 
of the 1st DUP. This process appears to be universal and inde- 
pendent of the stellar environment : it affects more than 95% of 
low-mass stars (Charbonnel & do Nascimento l 1 998> . whether 
they belong to the field, to open or globular clusters. Some indi- 
cations of such a process have also been detected in the bright- 
est RGB stars of external galaxies like the LMC (Smith et al. 

and Sculptor (Geisler et al. l2005t . Its signatures m terms 
of abundance anomalies are clear : the Li and '^C abundances 
as well as the '^C/'^C ratio drop while the '"^N and the ""O/'^O 
ratio increase. These data attest the presence of a non-standard 
mixing process connecting the stellar convective envelope with 
the external layers of the HBS where CN-burning occurs. Last 
but not least, they indicate that the effects of this process on sur- 
face abundances appear when the star reaches the RGB bump. 

Why should the RGB bump be such a special evolutionary 
point in the present context? After the completion of the 1st 
DUP, the CE retreats leaving a discontinuity of mean molecu- 
lar weight (or //-barrier) at the mass coordinate of its maximum 
penetration. Subsequently, when the HBS eventually crosses 
this discontinuity, the star suffers a structural re-adjustment due 
to the composition changes (more H is made available in the 
burning shell). The resulting alteration of the energetics causes 
a momentary decrease of the stellar luminosity. This results in 
a higher probability for finding a star in this brightness bin, and 
translates into a bump in the luminosity functions of globular 
clusters. Standard theory and observations nicely agree on the 
size and the location of the bump in the HRD (e.g., Zoccali et 
al. rT999ll . 

As for mixing in radiative stellar interiors, it was suggested that 
the discontinuity of molecular weight left by the 1st DUP could 
inhibit any extra-mixing between the base of the convective en- 
velope and the HBS. After the bump, the ju-gradients are much 
smoother in this region, permitting some extra-mixing to occur 
(Sweigart & Mengel [W7^ Charbonnel [T995l Charbonnel et al. 
[T998l l. 
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2. From the pioneering worl^ on stellar rotation to 
the present treatment of the transport of 
angular momentum and chemicals 

Sweigart & Mengel (119791 hereafter SM79) investigated the 
possibility that meridional circulation might lead to the mixing 
of CNO-processed material in RGB stars. Though the physics 
of rotation-induced mixing invoked at that time was very crude, 
this pioneering work has magnificently settled the basis of a 
complex problem. SM79 discussed in great details the problem 
of ju-gradients which were known to inhibit meridional circu- 
lation (Mestel fl 9531119571 . For mixing to be efficient, they un- 
derlined the necessity for the radiative zone separating the CE 
from the HBS not to present significant molecular weight gra- 
dients. The other crucial point made by SM79 concerned the 
importance of the angular momentum (hereafter AM) distribu- 
tion within the deep CE of RGB stars on the resulting mix- 
ing of CNO processed material. Indeed, beyond the 1st DUP, 
the angular velocity of a radiative layer near the HBS depends 
sensitively on how much AM has been deposited by the re- 
treating CE. SM79 investigated two extreme cases, namely a 
constant specific angular momentum and a uniform angular ve- 
locity within the CE. Substantial CNO processing of the enve- 
lope could be obtained with plausible MS angular velocity only 
when the inner part of the convection envelope was allowed to 
depart from solid body rotation. As we shall see in this paper, 
our lack of knowledge of the distribution of angular momentum 
within the CE of giant stars remains one of the weakest points 
of our global understanding of rotation-induced mixing at this 
phase. 

Rotational transport processes cannot be simply reduced to 
meridional circulation. Once established, this large scale circu- 
lation generates advection of AM, and thus favours the devel- 
opment of various hydrodynamical instabilities. Zahn J1992> 
proposed a description of the interaction between meridional 
circulation and shear turbulence, pushing forward the idea of 
chocking the meridional circulation by //-gradients. Following 
these developments but using a simplified version of Zahn's 
description, Charbonnel {1995} re-investigated the influence of 
such a process in RGB stars. She conjectured that the combina- 
tion of the disappearance of the mean molecular weight gradi- 
ent barrier after the bump and the increase of mixing in the HBS 
has the proper time dependence to account for the observed be- 
haviour of carbon isotopic ratios and for the Li abundances in 
Population 11 low-mass giants. 

In these exploratory computations however, the diffusion 
coefficient for chemicals was derived from an assumed constant 
rotation velocity (in the radiative zone) on the RGB, and the 
transport of AM by hydrodynamical processes was not consid- 
ered. This is however of utmost importance in the understand- 
ing of the rotation-induced mixing (see the review by Maeder 
& Mevnet l2000t . 

Then, Denissenkov & Tout J2000I I applied the formalism 
of Maeder & Zahn J1998> to a typical globular cluster RGB 
star. However, this was done in a post-processing approach and 
thus, did not take into account the feedback of mixing on the 
stellar structure. Considering the angular velocity at the base 
of the CE as a free adjustable parameter and treating the trans- 



port of chemicals only beyond the bwnp as in the aforemen- 
tioned works, they obtained large diffusion coefficients able 
to reproduce not only the Li, C and N abundance anomalies 
at the bump, but also the O-Na and Mg-Al anti-correlations. 
These two features being clearly primordial, Denissenkov & 
VandenBerg J2003t revised these results. They simplified their 
previous approach, letting AM evolve only due to structural 
readjustments (no rotational transport), and derived a diffusion 
coefficient to be applied to the chemicals beyond the bump. 
They considered the obtained mixing rate to have ''the correct 
order of magnitude" , even though it is too low by a factor of 7 
to reproduce the observational data for Pop I stars. 

In the present paper we propose a self-consistent approach 
of rotational-mixing in low-mass RGB stars. We define here as 
self-consistent a model in which the transport of angular mo- 
mentum and of chemicals is coupled to the evolution of the 
star from the zero age main sequence on. If meridional circu- 
lation and shear-induced turbulence are the only transport pro- 
cesses of AM considered, the assumptions made on the rota- 
tion profile concern solely the initial condition, i.e. the rotation 
profile at the ZAMS (assumed to be uniform) and the rotation 
regime in the convective envelope. At each evolutionary step 
we thus compute the new rotation profile together with the as- 
sociated transport coefficients resulting from structural read- 
justments and transport processes associated with rotation. The 
abundance profile of each chemical is then modified under the 
effect of both mixing and nuclear reactions. In such a proce- 
dure, the stellar structure "reacts" to rotational-mixing. 
We discuss the effects of rotation in RGB stars taking into ac- 
count the latest prescriptions for the transport of AM in stellar 
interiors, and the associated mixing processes. We describe the 
physical inputs of our models in §|3land §|3 and their effects on 
the angular velocity profiles and diffusion coefficients in §|5] In 
§|6lwe present the results for structural evolution and surface 
abundance variations of rotating low-mass Pop II stars from the 
ZAMS to the upper RGB. We then give a summary of our main 
results, and discuss them in relation with previous works in §0 
and propose new investigation paths in §|8l 



3. Physical inputs for the evolution of rotating 
stars 

3. 1 . Standard input physics 

The models presented here were computed with STAREVOL 
V2.30, and the reader is referred to Siess et al. (I2000I I and 
Palacios et al. ( |2003, paper I^ for a detailed description. Let us 
recall the main inputs. 

The nuclear reaction rates have been updated using the ver- 
sion 5.0 of the nuclear network generator NetGen available at 
lAA (http://astropc0.ulb.ac.be/Netgen). By default 
the adopted rates are : NACRE (Angulo et al. 19^ for charged 
particles, Bao et al. J2000> for neutron capture rates, Horiguchi 
et al. dl 9961 1 for experimental beta decay rates and Caughlan & 
Fowler J1988> otherwise. 
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For the radiative opacities, we use the OPAL tables^ above 
8000 K (Iglesias & Rogers I1996> and at lower temperatures 
the atomic and molecular opacities of Alexander & Fergusson 
J1994> . The conductive opacities are computed from a modified 
version of the Iben J1975> fits to the Hubbard & Lampe ( I1969> 
tables for non-relativistic electrons and from Itoh et al. ( I1983> 
and Mitake et al. Iil984t for relativistic electrons. 

The equation of state is described in detail in Siess et 
al. (2000) and accounts for the non ideal effects due to coulomb 
interactions and pressure ionization. The standard mixing 
length theory is used to model convection with ffMLT = 1 -75 
and the atmosphere is treated in the gray approximation and 
integrated up to an optical depth r - 5 x 10"^. 

3.2. Transport of angular momentum 



The evolution of AM and chemical species follow Zahn J1992> 
and Maeder & Zahn J1998> . Meridional circulation and turbu- 
lence induced by the secular shear instability are the two trans- 
port mechanisms considered here. Within this framework, the 
transport of AM obeys an advection/diffusion equation 



1 d , 4 . 15/4 5Q\ 



(1) 



where p, r and Q have their usual meaning, v,. is the vertical 
component of the turbulent viscosity associated with the shear 
instability. [/, is the vertical component of the meridional cir- 
culation velocity, which, assuming shellular rotation"*, is given 
by: 



M, 



■(En + Ef,) 



(2) 



where M« = M(l - j^^)- En and Zs^ depend respectively on 

the relative horizontal variation of the density, - p - 5^^' 

and of the mean molecular weight, ^ - j; - Detailed expressions 
for these quantities are given in Appendix A. The evolution of 
A depends on the competition between the vertical advection 
of a mean molecular weight gradient and its destruction 
by horizontal diffusion D/, (see § I3.4> 



dA dlnu 6 



dt 



dr 



This equation is obtained under the assumption that 



D, 



Di, 



(3) 



(4) 



Dy is the vertical turbulent diffusion coefficient (see § 13 .4> . 
and Ch are the characteristic distance scales in the vertical and 
horizontal directions respectively. As in paper I, the transport 



^ The OPAL tables used are based on a solar-scaled chemical mix- 
ture with possible enhancement of C and O, but do not include alpha- 
elements enhancement. 

"* The condition of shellular rotation is satisfied when turbulence is 
highly anisotropic and ensures CI x Q.(P). In that case, Zahn's for- 
malism may be applied strictly. Otherwise, it represents a first order 
approximation. 



of AM is computed by solving 5 first order differential equa- 
tions (Eq.|3]+ 4 equations resulting from the splitting of Eq.[0 
with a Newton-Raphson relaxation method. The upper bound- 
ary condition on A has however been modified compared to 
paper I, and we used the continuity equation at the base of the 
CE rather than assuming A = 0. 

Effects of yu-currents (Ef^ term in Eq. |2} are taken into ac- 
count consistently from the ZAMS up to the upper RGB. 

3.3. Transport of chemicals 

In the presence of strong anisotropic turbulence, Chaboyer & 
Zahn (1992) showed that the vertical advection of chemicals by 
a large scale circulation combined with strong horizontal diffu- 
sion produces a vertical effective diffusivity Deff (see § |^}. 
The vertical transport of a chemical species ; of concentration 
Ci can thus be described by a pure diffusion equation: 



dc, . i a r T n i a 
P— = Ci + — T- rpUipCi + — 
df ^ — „ — ' r~ or ^ J Or 



1 d 



1 d 



2 n dci 
r pDiot — 
or 



(5) 



nuclear 



atomic diffusion 



macroscopic processes 



Uip is the atomic diffusion velocity of the element with respect 
to protons, and Dtot is the total macroscopic diffusion coeffi- 
cient, and is the sum of the effective diffusion coefficient (|8} 
and of the vertical turbulent diffusion coefficient Dy (see § 13. 4> . 

The diffusion equation (|5jl is solved for each of the 53 
species considered in the code considering c, - (that is, diffu- 
sion and nucleosynthesis are decoupled). Here again, we used 
a Newton-Raphson method, as for the structure and AM trans- 
port equations. 



3.4. Diffusion and viscosity 

Let us briefly recall the various formulations used for the diffu- 
sion coefficients entering Eqs. Q, Q and (|5|i. 

• v,„D,, 

As in paper I, we assume that the secular shear instability 
dominates and that vertical shear eventually becomes turbu- 
lent in the radiative stellar interiors due to the low viscosity 
of the plasma (Zahn[1974). The development of turbulence 
is subject to the Reynolds criterion and sets in when 



Re = 



> Rec, 



Vmo\ + Viad 

where Vi-ad = '■^^ radiative viscosity, Vmoi is the 

molecular viscosity and Rec - 10 is the critical Reynolds 
number The shear instability obeys the Richardson crite- 
rion, and according to the classical formulation should set 
in when 

du 
dF ^ 

where Ric - 0.25 is the critical Richardson number. Here, 
we rather consider a modified Richardson criterion to take 
into account radiative losses and/or horizontal diffusion as 
described below. Several modified criteria have been pro- 
posed; in this paper we will compare results obtained with 
two of them. 
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- The first criterion assumes that thermal diffusion re- 
duces the stabilising effect of thermal stratification 
without affecting the chemical part. It leads to (cf. 
Maeder & Mey net 1996 1: 

8 Ric{rdQ.IArf - nJ; 



D, 



(6) 



where - Nj + NJ; is the Brunt-Vaisaala frequency 
and Kt is the thermal diffusivity. In this paper we will 
refer to this prescription as MM96. 
- The second criterion, which we used in paper I, also 
considers the erosion of the chemical stratification by 
the large horizontal diffusion (cf. Talon & Zahn lI99'7t : 



Dv -Vy--- 



(7) 



5Nl/(KT+Di,) + N'-/Dh' 
where D/, is the horizontal turbulent viscosity. In this 
paper we will refer to this prescription as TZ97. 
Let us further note that, due to the lack of a better prescrip- 
tion, we assume D,, = v,,. 

£>eff 

The effective diffusion coefficient solely appears in the 
transport equation of chemicals. In the approximation of 
highly anisotropic turbulence, it is the diffusive represen- 
tation of the effects of meridional circulation (Chaboyer & 
Zahn lI99'2t . and can be written as follows : 



£>eff = 



\ru,r 

30D,, 



(8) 



Vh, Dh 

The vertical turbulent diffusion coefficient D,, as well as the 
effective diffusion coefficient Dgfr depend on the horizontal 
component of the turbulent diffusivity D/,. No description 
of this diffusivity can be drawn from first principles, and 
its expression has to be parametrised. Assuming that the 
differential rotation on an isobar is small compared to unity, 
Zahn ( 1992 1 first proposed the following expression for v/,: 



Ch 



1 d(pp-U) UAXnr^Q. 



3pr dr 



2 dlnr 



= —\2V-aU\, (9) 



where C/, is a free parameter of order 1 which we used in 
our paper I, V = [p^'^U^ is the horizontal component 

of the meridional velocity and a - Meanwhile, 
some improvements have been achieved and more realistic 
prescriptions for the horizontal shear turbulent diffusivity 
including a dependence on the rotation rate were derived 
by Maeder ( t2003i and Mathis et al. (2004J : 



Vh 



^X" r [r Q(r) V'' [2V - at/])" , 



(10) 



where {x \ n ; k) - (40^ ; 3 ; 1) in Maeder's expression 
(with m = 1,3 or 5), and (1.610"*' ; 2 ; 0) in Mathis et 
al. J2004t . In Maeder J2003I I. the horizontal turbulent vis- 
cosity V/, is derived from a comparison between the dissipa- 
tion rate of turbulent energy by meridional circulation and 
the viscous dissipation rate, while in Mathis et al. J2004I I. it 
is derived from Couette-Taylor laboratory experiments (see 
Richard & Zahn [T^ . 

The horizontal turbulent diffusion coefficient D/, obtained 



by Eq. (I10> is larger than the one derived from expression 
Eq. (|9j, and is more consistent with the shellular rotation 
hypothesis. In the following, we will compare results ob- 
tained with these two prescriptions, which we refer to as 
Zahn92 and MPZOA (see TableQ. 

3.5. Rotation law in the convective envelope 



The formalism developed by Zahn J 19921 1 describes the trans- 
port of AM in radiative zones, and the rotation profile in the CE 
is defined by an upper boundary condition. Its choice is how- 
ever of prime importance since it determines the flux of angular 
momentum between these two regions. 

In the case of RGB stars, this will play an essential role as dis- 
cussed in §5.1. 

The interaction between rotation and convection is a long- 
standing and not fully understood problem, and despite the de- 



velopment of 3D numerical simulations (Ballot et al. 120041 
Browning et al. I2004> . the rotation profile within deep con- 
vective envelopes remains unknown. As already suggested by 
SM79, we may consider two limiting cases for the CE rotation 
law: 

1 . Uniform angular velocity ( solid body rotation ) 

This hypothesis has generally been assumed when mod- 



elling the evolution of rotating stars (Endal & Sofia 1 19781 
Talon et al. [T9971 Meynet & Maeder I2OOOI paper I). It 
is also the rotation regime obtained when describing AM 
transport in the CE by a diffusion equation using the dif- 
fusion coefficient derived from the MLT theory (Heger et 
al. 20001. Imposing solid body rotation in the CE is thus 
equivalent to assuming that the meridional currents are in- 
hibited in the presence of convection, and that the turbu- 
lent viscosity associated with convection is large enough to 
allow for instantaneous homogenisation of the angular ve- 
locity profile Q(r), as it is the case for chemicals (Endal & 
Sofia ll97'6t . This condition is also motivated by observa- 
tions of the solar convection zone (the sole star for which 
we have such informations at present), where radial differ- 
ential rotation is minute (Kosovichev et al. ll997> . 
2. Uniform specific angular momentum ( differential rotation) 



Already in the early 70's, Tayler ( I1973> addressed the ef- 
fects of rotation in stellar convective zones and came to the 
conclusion that meridional currents could develop and al- 
ter the rotation law in stellar convective regions. With cau- 
tion, he suggested that "it is possible that the asymptotic 
state [of rotation in a convective zone ] is closer to one 
of uniform angular momentum than uniform angular ve- 
locity". Some years later, Sweigart & Mengel J1979> pro- 
posed that rotation-induced mixing by meridional circu- 
lation could explain both the CNO abundance anomalies 
in RGB stars and the slow rotation rates observed in MS 
low-mass stars provided the radiative interior conserves 
its angular momentum during the first ascent of the giant 
branch, and the CE has constant (djcE/df = 0) and uni- 
form {djcE/dr = est) specific angular momentum. More 
recently, Denissenkov & Tout ( 120001 investigated both pos- 
sibiUties for rotation in the CE of a low-mass RGB star, and 
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Table 1. Rotation input physics of the models. D/, : MPZ04 
refers to Mathis et al. (2004) whereas Zahn92 refers to Zahn 
(1992). Dy : TZ97 refers to Talon & Zahn (1997) and MM96 
refers to Maeder & Meynet (1996). Uniform specific angular 
momentum in the CE was only applied beyond the turn-off. 
All models were computed with uniform angular velocity in 
the CE during the MS phase. Braking, if applied, begins at the 
ZAMS following a Kawaler (1988) law calibrated in order to 
lead to equatorial rotational velocity of v ^ 10 km.s"' when 
they reach the Hyades age (this occurs for Teff ~ 6300 K) early 
on the main sequence. This velocity is typical of cool stars in 
the Hyades. The highlighted row indicates our reference model 
Ml. 



Model 


CE 
rotation law 


km.s"' 


Braking 


D,, 


/i-euiTents 


D, 


MO 















Ml 


D(r) = est 


5 


no 


MPZ04 


yes 


TZ97 


M2 


j{r) = est 


5 


no 


MPZ04 


yes 


TZ97 


M3 


j{r) = est 


5 


no 


MPZ04 


no 


TZ97 


M4 


j{r) = est 


5 


no 


MPZ04 


yes 


MM96 


M5 


i{r) = est 


5 


no 


Zahn92 


yes 


TZ97 


M6 


j{r) = est 


110 


yes 


MPZ04 


yes 


TZ97 



concluded in favour of uniform and constant specific angu- 
lar momentum in the CE. Investigating the rotation rates on 
horizontal branch (HB) stars of globular clusters. Sills & 
Pinsonneault J2000t showed that slow rotation on MS stars 
and V sin ; up to 40 km.s ' on the HB indicate that a non- 
negligible amount of AM is preserved in the stellar interior 
between the MS turn-off and the HB. They proposed that 
this could be achieved assuming uniform specific angular 
momentum in the CE during the RGB phase, a conclusion 
again similar to that drawn by SM79 and more recently, by 
Chaname et al. (2004a, b). 

These studies indicate that the condition of differential ro- 
tation in the CE (i.e. jcB(r) - \r^O.{_r) - est) could be a key 
ingredient to derive a consistent history of the AM evolu- 
tion in low-mass stars, as well as to produce a high degree 
of mixing in the RGB interiors needed to explain part of the 
abundance anomahes observed at this phase. 



4. Numerical simulations 

We model a typical metal-poor globular cluster RGB star, with 
initial mass M - 0.85 Mq, initial helium content Kini = 0.248 
and [Fe/H] - -1.57 which corresponds to the metallicity 
of M13 stars (Sneden et al. I2004> . We take into account the 
a-enrichment expected at such metallicities, with [O/Fe] = 
[Ne/Fe] = [Mg/Fe] = -1-0.3 dex, and account for the odd-even 
effect on sodium, with [Na/Fe] = -0.3 dex. The initial ratio 
for the magnesium isotopes is -"^Mg Mg Mg = 87.5 : 
6 : 6.5, and is similar to the values determined by Yong et al. 
J2003> for the stars with the lowest ^^Mg and -''Mg in NGC 
6752, a globular cluster with [Fe/H]^ -1.6. The mass fractions 
of the other elements are in solar system proportion, and scaled 



for their sum to be equal to unity. This corresponds to a metal- 
licity of Z = 8 10 

Mass loss is included from the ZAMS on in all our models. 
We use the empirical Reimers J1975t formula with a metallicity 
scaling : 

M = -3.9810-'^'77^ J:^ Moyr-' (11) 
M ^ Z0 

with rj - 0.5. At each time step, the associated AM losses are 
taken into account and the adjusted total AM is conserved when 
applying rotational transport. 

All models are computed assuming uniform angular veloc- 
ity in the CE during the MS as indicated by the solar case. 
When considered, the hypothesis of uniform specific angular 
momentum in the CE is applied beyond the turn-off, i.e. when 
the convective envelope begins to deepen. We will compare two 
different rotational histories, namely models that were slow ro- 
tators already on the ZAMS and models with an initially larger 
velocity typical of ZAMS Pop I stars, but which experience 
magnetic braking in their early evolution, as expected from 
solar-type stars. This braking results in similar surface veloc- 
ities at half-way on the main sequence (when Xc x 0.5) for ini- 
tially fast and slow rotators. Considering their important role 
in shaping the rotation profile (see paper I), the so-called /i- 
currents (E^ term in Eq. |5J are taken into account in all our 
rotating models except M3. In all cases the initial rotation pro- 
file at the ZAMS is defined by Q(r) = vm/R* - est V r. 

We finally underline that Eq. Q is solved in all its complex- 
ity for all the rotating models presented here from the ZAMS 
up to the upper RGB. 

Tabled lists the characteristics of the models that we have 
computed using different prescriptions for the input physics. 
For our reference model (Ml), we consider the following set 
of parameters/physical ingredients: uniform angular velocity 
in the convective regions at all times, initial surface velocity 
of 5 km.s ' on the ZAMS and no braking applied, Mathis et al. 
( l2004l : Eq.[lO|l and Talon & Zahn ( fT997l : Eq.Ell prescriptions 
for the horizontal v/, and vertical v,, turbulent diffusion coeffi- 
cients respectively. 

Figure IT31 presents the Hertzsprung-Russell diagram and 
the evolution of the surface luminosity and temperature for the 
models computed. It will be discussed in more details in § 16.11 
The evolutionary points on which we will focus in the follow- 
ing sections are marked on the evolutionary path and corre- 
spond to the turn-off (A), the end of the first DUP (B), the bump 
(C) andL = 3IOL0 (D). 

5. Testing the physics of rotation 

In this section, we analyse the impact of the different physical 
inputs included in the models presented in Table[n 

5.1. Rotation in tine convective envelope 

In § 13.51 we mentioned the uncertainty regarding the rota- 
tion law in the convective envelope of a giant star Previous 
studies also provide some hints that the evolution of the AM 
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Fig. 1. Upper panel: Hertzsprung-Russell diagram for the different test models presented. All sequences start with Mzams — 
0.85 Mq, [Fe/H] = -1.57. Evolutionary tracks for models Ml and M3 are superimposed. Open symbols labelled A, B, C and D 
correspond to the turn-off, the end of the 1 st DUP, the bump and L - 3 1 L© respectively. Lower panels: Evolution of luminosity 
(left) and effective temperature (right) as a function of time beyond the turn-off. 



(and chemicals) distribution within the radiative interior might 
strongly depend on the rotation regime at the base of the con- 
vective envelope. This aspect can be studied by comparing 
models Ml and M2, which solely differ by the applied rotation 
law in the CE beyond the turn-off. 

Profiles of the angular velocity Q., of the horizontal density 
fluctuations O { - % - \^dQ.ldr) and of the specific angular 
momentum j at four different evolutionary stages for models 
Ml, M2 and M6 are presented in Figs. |2] and |3l In Fig.|2j as 
well as in several other figures in this paper, quantities are plot- 
ted against 6M instead of Mr- 6M is a relative mass coordinate 
allowing for a blow-up of the radiative region above the HBS, 
and is defined as 

Mr - Mhbs 



6M = 



Ml 



BCE 



-M„ 



(12) 



6M is equal to 1 at the base of the convective envelope and 
at the base of the HBS (where X = 10"^). Typically in our 
models, the nuclear reactions occur between 6M = 0.2 and 6M 
- 0.5, a large mean molecular gradient being associated with 
the layers of maximum energy production at 6M - 0.2 

By assumption, the rotational evolution is identical in mod- 
els Ml and M2 up to the turn-off. Beyond, the evolution of 



the AM distribution is dominated by structural readjustments 
of the star becoming a giant. The degree of differential rotation 
in the radiative zone globally increases with time, leading to a 
rapidly rotating core and a slowly rotating surface (Q rises by 
3 to 4 orders of magnitude between the base of the CE and the 
edge of the degenerate He core in both models). 
At the turn-off, the base of the CE rotates at the same veloc- 
ity in models Ml and M2. As the star crosses the Hertzsprung 
gap and approaches the Hayashi line on its way to the red giant 
branch, the core contracts while the outer layers expand and 
are thus efficiently slowed down. For model Ml, the velocity 
of the inner shells of the envelope also decreases substantially 
because of the solid body rotation of the CE (dQc£/dr = 0), 
and the AM attached to the convective envelope concentrates 
in the outer layers. At the end of the 1st DUP the angular ve- 
locity at the base of the CE, Hbce, decreased by a factor of 
65. When the CE withdraws in mass, it leaves behind radiative 
shells with low specific angular momentum that rotate slowly. 
The resulting differential rotation rate is low in this region as 
shown in the and Q. profiles. 

For model M2, during the 1st DUP the surface layers slow 
down but the uniformity of the specific angular momentum pre- 
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Fig. 2. Angular velocity Q and horizontal density fluctuations 9 oc dD./dr inside models Ml (solid lines), M2 (dotted lines) and 
M6 (dashed lines) at the four evolutionary points top-labelled on the graph and reported on Fig l3.5l £1 and 9 are plotted against the 
scaled mass coordinate 6M which allows a blow-up of the region of interest {6M = at the base of the HBS, and 6M = 1 at the 
base of the CE). The long dashed lines on the upper pannels represent the D. profiles that one gets when the angular momentum 
evolution in model M2 only results from the structural readjustments beyond the turn-off (i.e. no AM transport). 



vents the shells of the inner CE from decelerating abruptly and 
ensures the concentration of AM in this denser region. At the 
end of the 1st DUP, Q.bce has decreased only by a factor of 3. 
As the CE withdraws in mass, its deeper shells retaining the 
larger part of the CE angular momentum fall back into the un- 
derlying radiative zone, where they conserve their large angular 
velocity and strong differential rotation. 



Figure|3shows that after the completion of the 1st DUP, the 
assumption of uniform specific angular momentum within the 
CE also translates into a quasi-constancy of the level of spe- 
cific angular momentum in this region (dotted and dashed lines 
in panels (B), (C) and (D)). This is a natural expectation since 
the convective envelope of a red giant represents more than 80 
% of the stellar radius and retains most of the angular momen- 
tum. The assumption of constancy of jcE under a regime of 
uniform specific angular momentum remains however an ap- 
proximation since transfer of AM occurs between the interior 
and the envelope. Our models indicate that the specific angular 
momentum in the envelope indeed slightly increases with time. 



The difference of the rotation profiles encountered in mod- 
els Ml and M2 below the CE affects the diffusion coefficients, 
in particular D,,, which scales as (rdQ/dr)^. 

Figures 0] and |5] present the profiles of the diffusion coef- 
ficients entering the transport equations of AM and chemicals 
(Eqs.^andlS} at the end of the 1st DUP, at the bump and at 
L = 310 Lq for models Ml and M2, respectively. One may 
notice that Z);, is everywhere much larger than £),, (see also 
Figs. |6j|^ and This validates the shellular rotation hypothe- 
sis. 

- Model Ml 

At the end of the 1st DUP (left panel in Fig.|3|i, differen- 
tial rotation is negligible below the CE in the region be- 
tween 6M ^ 0.5 and 6M - 0.91. There, the vertical tur- 
bulent viscosity is smaller than its critical value given by 
i^i'.crit - ^^cViad - 3000, and turbulence does not develop. 
The low rate of differential rotation quenches the turbulent 
transport and creates a gap that disconnects the CE from 
the regions where nucleosynthesis occurs. Below 6M ^ 0.5 
the differential rotation is larger and the shear turbulence 
dominates the transport down to the region {6M ^ 0.1) 
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Fig. 3. Profile of the logarithm of the specific angular momentum log(7) inside models Ml (solid lines), M2 (dotted lines) and 
M6 (dashed lines) at the four evolutionary points top-labelled on the graph and reported on Fig. 13. 51 as a function of the internal 
mass coordinate Mr in units of M©. 



where both meridional circulation and shear turbulence 
are very efficiently hindered by the high mean molecular 
weight barrier associated with the hydrogen burning (see 
upper panels of Fig. llOt . At the bump and beyond, shear- 
induced turbulence does not develop and the effective dif- 
fusion coefficient remains negligible below 6M ^ 0.85, so 
that no modification of the surface abundance pattern by ro- 
tational mixing should be expected beyond this point (see 

§0- 

- Model M2 

When assuming uniform specific angular momentum in 
the CE, the picture is very different (Fig. O. The larger 
differential rotation rate allows the shear flow to become 
turbulent almost everywhere between the base of the CE 
and the HBS. Shear turbulence dominates the transport 
of chemicals across the entire radiative zone, the effec- 
tive diffusion coefficient associated with meridional circu- 
lation being always much smaller in the whole radiative re- 
gion. Anticipating the results of § 16.21 we can already see 
from Fig. ^1 that Aot does not rise above lO^cm^.s"' in 
the outer HBS (around 6M = 0.2-0.3) even at the bump, 
which is much lower than the parametric "canonical mix- 



ing rate" of 4.10^cm-.s ' deduced from observational con- 
straints (Denissenkov & VandenBerg |2003t . 



The profile of Dgff in Fig |5l appears to be quite ragged be- 
low 6M ^ 0.5, in particular on panel (D). Each of the bump in 
this region is associated with an inversion of the meridional cir- 
culation velocity, which can be positive or negative. Although 
the meridional circulation can actually present various cells, 
these particular features are due to numerical instabilities oc- 
curring in the regions where the mean molecular weight gra- 
dients are non-negligible. When the effects of //-currents are 
included (terms depending on // or A in Eq. |2j, the numerical 
system is highly non-linear and may be sensitive to numerical 
parameters such as the spatial and temporal resolutions. 
The global amplitude of Dgg (and Dy) in these regions of large 
;7ZM-gradients, together with the decrease in the diffusion co- 
efficients in the nuclearly active shells are however robust re- 
sults. Just below the CE, the first bump in the Detf profile is also 
a robust feature, associated with the Gratton-Opik meridional 
circulation cell. These remarks apply to all our models. 
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Fig. 4. Profiles of the different diffusion coefficients at the end 
of the first dredge-up (left panel) and at the bump (right panel) 
for our reference model (Ml in Table 1). 




Fig. 5. Same as Fig.l^for M2 in Table 1. M2 only differs from 
Ml by the adopted differential rotation law in the convective 
envelope 

5.2. Impact of the ZAMS rotation velocity 

We investigated the impact of the initial rotational velocity 
considering an originally slow (model M2) and a fast rotator 
(model M6) on the ZAMS, without changing any of the other 
physical parameters. We have applied to model M6 the same 



treatment as to Pop 1 stars (see paper I), namely a strong mag- 
netic braking according to the Kawaler J1988> prescription so 
as to get a surface equatorial velocity lower than 10 km.s"' 
before the central hydrogen mass fractions gets lower than 0.5. 
At the turn-off the surface velocities of models M6 and M2 are 
thus quite similar (see Table|3}. 

The point under scrutiny in this section is to determine whether 
this strong braking, which triggers strong turbulence in the ra- 
diative zone during the main sequence, has an incidence on 
the angular momentum distribution and on the diffusion coeffi- 
cients beyond the turn-off. 

In the present study, we did not investigate the case of fast rota- 
tors at the turn-off, this configuration being ruled out by the ve- 
locities available from observations of low-mass turn-off stars 
in globular clusters (Lucatello & Gratton E003t . 

Figure |2| presents the degree of differential rotation {9 
profiles) and the angular velocity inside models M2 and M6 
(dotted and dashed lines respectively) at different evolutionary 
points. Model M6 rotates globally faster than model M2 at the 
turn-off, and this difference is maintained during the evolution. 
Indeed, although model M6 undergoes a very efficient braking 
on the early MS, it is only braked down to 5.7 km.s ' at the 
turn-off, compared to 3.88 km.s ' in model M2. In spite of 
these different surface rotation rates, the profiles of Q and 
are quite similar in the radiative zone of these models during 
the RGB phase (Fig. |3}. Consequently, the turbulent diffusion 
coefficients (see Eq. are not very different during the RGB 
phase, as can be seen on Figs.|5]|6land [TO] 



This comparison shows that the structural readjustments 
(induced the 1st DUP) efficiently redistribute the AM through- 
out the star beyond the turn-off. The resulting rotation profile 
after the completion of the 1st DUP is almost solely determined 
by the star's total angular momentum at the turn-off and by the 
CE rotation law. 



5.3. /j-currents and horizontal turbulence 

The mean molecular weight affects the transport of chemicals 
via its gradient V,, and its relative horizontal variation A (see 
Eas.l2land lA.3> . /i-gradients inhibit the efficiency of meridional 
circulation, and /i-currents are inhibited by strong horizontal 
turbulence. In paper 1 we emphasised the importance of these 
terms in establishing the differential rotation profile during the 
MS for Pop 1 stars. 

In order to better assess the role of /i-currents in Pop 11 stars, 
and to determine their sensitivity to the horizontal turbulence 
description, we have computed a series of models for which 
we alternately use the Zahn92 and the MPZQA prescriptions for 
Dh, assuming both Ef^ - and Efi + 0, and vzams - 5 km.s"' 
and Vzams = 110 km.s"' for each case. The properties of these 
models are summarised in Table|2] 
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Fig. 7. Profiles of the vertical component of the meridional circulation velocity U,-, angular velocity Q and total diffusion coef- 
ficient Dtot inside 0.85 M©, [Fe/H] = -1.57 models at three stages of their evolution (as indicated by their central H content Xc : 
= 0.6 <=> f ~ 2 Gyr, Xe = 0.15 <=> f 8 Gyr and Xe - 0.01 <=> f ~ 10 Gyr). First 3 columns on the left Profiles for models that 
are slow rotators on the ZAMS {vzams - 5 km.s ') and undergo no braking. Loif 3 columns on the right Profiles of initially fast 
rotating models on the ZAMS {vzams = 110 km.s"') and undergoing strong magnetic braking during their early MS evolution. 
Here, sohd and dashed lines are superimposed and indistinguishable from one another. Details on input parameters are given in 
Table 2. 




Fig. 6. Same as Fig. |3 for M6. In this model, vzams = 
110 km.s ' and after efficient braking through according to a 
Kawaler law , uto = 6.4 km.s"' 

. Uniform specific angular momentum was assumed in the CE. 

Main sequence 

Let us first comment on the MS evolution. The profiles of Q., 
Ur and Dtot are displayed on Figure Q at three different times 
on the MS fsee Tahlel2lfor detailed descrintion of each curve"). 



Table 2. Input parameters for the models discussed in ^5.31 
Some of the models appearing in this table are also listed in 
Tab. 1 . Models Ma, Mb, Mc, Md and Mh have only been com- 
puted up to the turn-off. In all models, solid-body rotation in 
the CE has been assumed, and D,, is given by TZ97. 



Model 


i'ZAMS 


Braking 


Di, 


/i-cuiTents 


line style 












in Fig. 7 


Ma 


110 


yes 


MPZ04 


no 


solid 


Mb 


110 


yes 


MPZ04 


yes 


dotted 


Mc 


110 


yes 


Zahn92 


yes 


shoit-dashed 


Md 


110 


yes 


Zahn92 


no 


long-dashed 


M3 


5 


no 


MPZ04 


no 


solid 


M2 


5 


no 


MPZ04 


yes 


dotted 


M5 


5 


no 


Zahn92 


yes 


shoit-dashed 


Mh 


5 


no 


Zahn92 


yes 


long-dashed 



Slow rotators on the ZAMS (Fig. 7, left panels). 

Their slow rotation together with the negligible structural 
readjustments occurring during the MS does not favour any 
mechanism able to trigger steep rotation profiles. As a result, 
the transport of both angular momentum and chemicals by 
rotation-induced processes is inefficient. This is independent of 
the choice for the D/, prescription and of the introduction or not 
of the //-currents. The Q-profile as well as the total diffusion 
coefficient for chemicals is similar in all cases, yu-gradients 
in the radiative interior are small, which makes the effects of 
//-currents negligible during the MS. 



12 



Palacios et al.: Rotational mixing in low-mass stars. 11 



Fast rotators on the ZAMS (Fig. 7, right panels). 

Turning now to the fast rotators on the ZAMS, the strong 
braking at the beginning of the MS allows the build-up of 
steeper Q-gradients, and shear-turbulence can develop across 
the radiative zone akeady during the main sequence. This leads 
to efficient transport of both angular momentum and chemicals, 
similarly to what is obtained in their Pop 1 counterparts. The //- 
gradients are larger in this case than for slow rotators, and we 
can expect /i-currents to have the same effects on the rotation 
profile as those observed in the Pop 1 stars. It is actually the 
case when considering the Zahn92 prescription for D/, which 
was used in paper 1. As shown on Fig when the /^-currents 
are taken into account (E^ + 0; short-dashed lines), the de- 
gree of differential rotation reached near the end of the MS is 
substantially larger than when these terms are neglected (long- 
dashed lines). The yu-currents also limit the extent of the shear 
unstable region at the centre, whereas when /s^ = 0, turbulence 
is free to develop across the entire radiative interior (third col- 
umn, Fig-EJ- 

Using the MPZ04 for D/, leads to a different conclusion. 
This prescription produces much larger values for the horizon- 
tal turbulence D/, than the former one. As the /^-currents {Ef, 
term in Eq.|2l see also EQ. IA.3l in Appendix) are generated by 
the horizontal variation of the mean molecular weight, they are 
much reduced by this choice. In the case we are considering, 
they are reduced to an insignificant level, i.e. they remain small 
compared to the Q-currents (term E^ in Eq.|2 see also Eq. lA.U 
and do not affect the rotation profile nor the meridional circu- 
lation velocity. 

In short, on the MS, the yU-currents have negligible effects 
on Pop 11 low-mass stars that are slow rotators, independently 
of the prescription used for D/,. This is exclusively due to the 
slow rotation. On the other hand, the ju-currents may affect the 
building of the rotation profile during the MS evolution of low- 
mass Pop 11 stars undergoing strong magnetic braking if the 
horizontal turbulence is not too large to prevent any significant 
variation of the mean molecular weight to develop. 

Red Giant Branch 

Comparing model M2 with model M3 and M5 provides respec- 
tive clues on the effects of the yU-currents (M2 versus M3) and 
the horizontal turbulence description (M2 versus M5) on the 
transport of chemicals beyond the turn-off. 
Figure^lpresents the total diffusion coefficients for the chem- 
icals in all our rotating models at evolutionary points (B), (C) 
and (D). Dtot in models M2, M3 and M5 (dotted, long-dashed 
dotted and long-dashed lines respectively) is essentially the 
same at the different evolutionary points presented. This indi- 
cates that /i-currents do not affect the transport of chemicals 
on the RGB. On the other hand, the choice of the prescrip- 
tion for horizontal turbulence has small effects on Dtot. This is 
due to the fact that £)/, essentially affects the effective coeffi- 
cient Deff, which is much smaller than D,, after the 1st DUP in 
case Zahn92 is used, and remains small for MPZ04. Let us add 
that in general, Dh,MPZ04 > Dh,zahn92, and that for Dh,zahn92, the 




0.5 1 0.5 1 0.5 1 
6M 6M 5M. 

Fig. 8. Same as Fig.|3]for M5. This model is similar to M2 but 
was computed using the Zahn (1992) expression for D/,. 



shellular rotation hypothesis (i.e. D/, » D, ) is violated in some 
regions (see panel (D) Fig.|S}. 

5.4. Vertical turbulent diffusion coefficient 

In order to compare our results with the calculations of 
Denissenkov & Tout (I2000> . we also computed a model (M4) 
using the MM96 prescription for v,, (Eq. |6jl. The main differ- 
ences between this criterion and the one derived by Talon & 
Zahn J1997> were already exposed in § 13.41 In Fig.|9]we present 
the diffusion coefficients for model M4 at the end of the 1st 
DUP at the bump and at L ^ 310 Lq. Contrary to models M2, 
M3, M5 and M6, for which we also assumed uniform specific 
angular momentum in the CE beyond the turn-off, by the end 
of the 1st DUP, the shear instability has not developed between 
the base of the CE and the edge of the HBS in model M4. 
Indeed Eq. (|6j shows that the shear must be larger than NJ^ for 
the instability to develop. In the case of giants, it is only after 
completion of the 1st DUP, where the retreating CE leaves a 
chemically homogeneous region (A^^ = 0), that the shear insta- 
bility may set in. However the CE remains disconnected from 
the nucleosynthesis regions until the large jU-barrier left by the 
DUP is erased at the bump. This justifies the assumption of 
Denissenkov & Tout J2000> . that considers rotational transport 
of chemicals only from the bump on. 

However, Meynet & Maeder ( I1997> showed that using 
this strong criterion prevents mixing from occuring in mas- 
sive, fast rotating stars, in contradiction with observational ev- 
idence. The same conclusion that was also reached by Talon et 
al. J1997> . and motivated the authors to introduce the erosion 
of the ju-gradient by horizontal turbulence. For the same rea- 
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Fig. 9. Same as Fig. |4]for M4. This model is similar to M2 
except that the prescription for the vertical turbulent diffusion 
coefficient is that of Maeder & Meynet (1996). 



son, Maeder (I1997> also developed a modified shear criterion 
to reduce the efficiency of mean molecular weight gradients. 

Being in the same framework (i.e. study of transport asso- 
ciated with meridional circulation and shear turbulence), and 
in absence of strong observational evidence, it is not justified 
to change these prescriptions for the particular case of low- 
mass RGB stars. It was SM79 who first suggested that prior 
to the bump, mixing should be hindered by the //-barrier left at 
the end of the 1st DUP, and that surface abundance variations 
should not be expected before this evolutionary point. We will 
show however in the following that a mixing process eroding 
the ju-gradient does not necessarily alter the surface abundance 
pattern prior to the bump. 

5.5. Dominant process for the transport of angular 
momentum and chemicals 

In the framework of this paper, we consider two transport pro- 
cesses, namely meridional circulation and turbulence induced 
by the secular shear instability. Whereas shear-induced turbu- 
lence is always described as a diffusive process, meridional cir- 
culation appears as an advective process for AM and as a diffu- 
sive process for the chemical species (Eq.^^nd Eq.lsjrespec- 
tively). The relative importance of these processes can thus be 
different depending on whether we consider AM or elements 
transport. 

In the case of chemical species. Figs. |?}|6] and indicate 
clearly that shear-induced turbulence is the dominant transport 
process beyond the completion of the 1st DUP when uniform 
specific angular momentum is assumed in the envelope. The 



Fig. 10. Lower panels Total diffusion coefficient for several 
models (line encoding is given in the figure) as a function of the 
reduced mass coordinate 6M at the end of the first dredge-up 
and at the bump. Upper panels Hydrogen mass fraction profiles 
in the same models. 



exception is model M5, where due to the smaller efficiency of 
the horizontal shear turbulence (D/, = Df,,zahn92), meridional cir- 
culation is large and still dominates the transport of chemicals 
at the end of the 1st DUP (Fig. E}. Higher on the RGB how- 
ever, turbulence recovers ffie upper hand in this model too. For 
model Ml, shear-induced turbulence can not be triggered in ffie 
radiative zone after the 1st DUP, and meridional circulation ac- 
counts for transport for chemicals. The value of Dgif remains 
however always smaller that the molecular viscosity v„,oi. 

In order to estimate the efficiency of the two processes 
responsible for the AM transport, we compare their relative 
characteristic timescales. The characteristic timescales for 
AM transport by meridional circulation and by shear-induced 
turbulence over a distance Ar can be estimated as T[/ ~ ^ and 

Ty ~ respectively. 

Tu « Ty at the end of the 1st DUP for all models but model 
Ml. This situation is subsequently modified as shells with a 
large angular momentum are incorporated from the inner CE 
into the radiative interior From the bump on, the efficiency of 
turbulence becomes thus similar to that of meridional circula- 
tion, with Ty ^ Tu ~ 10^ yr. For model Ml, r,, ^ Tf/ ^ 10^ yr 
from 6M ~ 0.1 to 0.8 at all times (evolutionary points (B), (C) 
and (D)). In the region just below the CE, for 6M 6 [0.8; 1], 
the degree of differential rotation is very small (see Fig. O so 
that the timescale for AM transport by turbulence is of the 
order of 10 Gyr, and meridional circulation dominates. 
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Between the turn-off and the completion of the 1st DUP, 
meridional circulation dominates the transport of AM in all 
our models, while the chemical species are mainly transported 
through turbulence. Beyond the 1st DUP, chemicals are trans- 
ported via shear-induced turbulence while the AM essentially 
evolves due to structural readjustments. Indeed at this phase, 
the advective (meridional circulation) and the diffusive (shear- 
induced turbulence) processes almost compensate each other, 
letting the Lagrangian term control the evolution of the angular 
velocity profile. In those models where turbulence can not de- 
velop, meridional circulation dominates the transport of chem- 
icals and AM (model Ml). 

The dominant transport process for AM can thus differ from 
the process controlling the transport of chemicals. It can also 
vary along the evolution. 

6. Signatures of mixing 

6.1. Structure evolution 

Through the transport of chemicals, and in particular those con- 
tributing to the nuclear energy production and the opacity, ro- 
tation may indirectly affect the evolution of stars in terms of 
lifetimes, luminosities and effective temperatures. Centrifugal 
forces can also affect the structure but this effect is negligible 
for the low rotation rates considered here. 

6.1 .1 . Main Sequence 

Figure l33] presents the Hertzsprung-Russell (HR) diagram for 
the models listed in Tabled Models MO to M3 can hardly be 
distinguished on the figure. Due to their small initial rotation 
velocity and the absence of braking to pump AM, models Ml 
to M3 present very weak differential rotation resulting in inef- 
ficient mixing during the MS phase (Aot ~ 10" - lO-' cm-s"'); 
their evolution as well as their chemical structure is thus only 
scarcely affected by rotation-induced mixing. Table |3]presents 
the main evolutionary characteristics of our models. The lumi- 
nosity and eflFective temperature at the turn-off, as well as the 
time spent on the MS are similar in slowly rotating (M1-M5) 
and standard (MO) models. 

Model M6 deviates from the standard and slow-rotating 
tracks. In this model, the strong braking applied during the 
first million years spent on the MS creates a large differen- 
tial rotation inside the star leading to an efficient transport of 
the chemicals (see Eq.0. In this case, as helium diffuses out- 
wards, the opacity is globally lower and the star consequently 
bluer and more luminous. As a result of fuel replenishment due 
to efficient rotational mixing, the H burning phase also lasts 
for a: 400 Myr longer in model M6, which is thus older at the 
turn-off compared to the other rotating models (Table |3}- 

For all our rotating models, the surface velocity at the end 
of the MS is lower than 6 km.s"' , in fair agreement with the up- 
per limits derived for globular clusters MS stars (Lucatello & 
Gratton l2003> . During the MS evolution, the surface rotation 
velocity remains almost constant in models Ml to M5. They 
are slowed from 5 km.s"' on the ZAMS to a; 4 km.s"' at 
the turn-off mainly due to the structural readjustment that be- 



Table 3. Main evolutionary features of the models presented 
in Tabled TO index stands for turn-off values. Mbce.dup is 
the mass coordinate associated with the deepest extension of 
the convective envelope during the first dredge-up. Afbump is 
the time needed for the entire HBS to pass through the fi- 
discontinuity left by the 1 st DUP. 



Model 





1 


2 


3 


4 


5 


6 


/to (Gyr) 


11.11 


11.20 


11.21 


11.21 


11.21 


11.21 


11.61 


L^o (Lq) 


3.35 


3.34 


3.35 


3.35 


3.35 


3.35 


3.55 


T'eff.TO (K) 


6534 


6547 


6543 


6543 


6543 


6543 


6582 


U^o (km.s-') 





3.88 


3.90 


3.89 


3.85 


3.95 


5.71 


A/BCE.DUP ( Mo) 


0.306 


0.303 


0.304 


0.306 


0.304 


0.307 


0.306 


^bump ( L©) 


107 


101 


109 


113 


103 


110 


114 




4899 


4919 


4897 


4887 


4910 


4896 


4897 


Arbump (Myr) 


4.87 


7.14 


7.75 


6.28 


5.75 


6.68 


7.73 



come important at the end of this phase. Model M6 undergoes 
magnetic braking on the MS so that its rotational velocity has 
already dropped below = 6 km.s ' when the model reaches the 
middle of the main sequence (i.e. for Xc - 0.5). This veloc- 
ity further decreases due to the efficient transport of AM in the 
radiative interior, and reaches 5.7 km.s ' at the turn-off. 

6.1.2. Red Giant Brancli 

The large increase in radius accompanying the deepening of 
the convective envelope during the 1st DUP, combined with 
global conservation of AM, leads to very efficient braking 
of the surface layers and spin-up of the core. All our rotat- 
ing models have surface velocities lower than 1 km.s"' at the 
end of the dredge-up (including model M6 for which we have 
vzAMS = 110 km.s '). In models Ml, M2, M4 and M6, for 
which the total diffusion coefficient for chemicals is larger than 
the molecular viscosity at the base of the envelope, the 1st DUP 
is deeper compared to the standard case. These variations of the 
depth of the 1st DUP (Tab.|3 remain however small and will 
not affect significantly the surface abundance patterns at this 
phase (see § 16.21 Fig. l6.2> . 

In the following we will refer to the bump luminosity as the 
luminosity of the model when the mass coordinate of the maxi- 
mum energy production inside the HBS is equal to Mbce.dup- 

Figures fin 1 121 and 1 131 present selected abundance profiles 
and the mean molecular weight gradient V^u - d In ju/dr, be- 
fore and after the bump inside models MO, Ml and M6 re- 
spectively. The peak amplitude ^ in the /^-gradient profile at 
Mr - Mbce.dup (indicated by an arrow on the left panels) de- 
creases with increasing degree of mixing. In the non-rotating 
model, this peak is a signature of the deepest penetration of the 
convective envelope during the first dredge-up, and the corre- 
sponding V^c value is ^ 10"" (Fig. II 11 1. In the rotating mod- 
els, the peak is spread out due to the ongoing transport of the 
chemicals and varies between V^ ^. ^ 3 ■ 10"'~ (Fig. I12li and 
V^.f ^ 5 ■ 10"'^ (Fig. [nil. The amplitude of V^.^ directly re- 
ffects the strength of the diffusion coefficient in this area (see 
Figs. |3and|6|l. According to Charbonnel et al. J 19981 1 the re- 
gions where > 1.5 x 10"'^ should not be affected by the 
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Fig. 11. Logarithm of the mass fraction profiles of H, '^C, '^C, 
''^N, '^O, ^^Na and of the mean molecular weight gradient = 
dlnfi/dr (bold line) as a function of the reduced mass inside a 
0.85 Mq, Z = 0.0005 standard stellar model MO. The panels 
represent the profiles of these elements in the HBS just before 
(i.e. before the HBS contacts the yu-discontinuity left by the 1st 
DUP) and after (i.e. once the mass coordinate at the base of the 
HBS is larger than Mbce.dup) the bump. In the left panel, the 
down-going arrow indicates the mass coordinate of maximum 
penetration of the CE during the 1st DUP. 

extra-mixing acting below the CE. Therefore, even the low 
value of c found in model M6 can prevent mixing to act 
freely below 6M ^ 0.5. 

Let us finally emphasise that although the amplitude of the dif- 
fusion coefficients may be locally increased due to numerical 
instabilities, and artificially lower the mean molecular gradient 
in that region, the reproducibility and constancy of the ju-barrier 
spread over indicates that this feature is not a numerical artifact. 

In order to evaluate the "observational" impact of rotational 
mixing on the luminosity function's bump, we have computed 
theoretical luminosity functions (hereafter LF) for each of our 
models (see Fig. ll4l i. These theoretical LF represent the time 
spent in each bin of magnitude V by models on the RGB with 
magnitudes lower than V = 16.6. We have chosen this cut-off" 
in order to be able to compare the models predictions with the 
observed luminosity function of the globular cluster Ml 3 as 
given by lCho et al. (2005)| 

In order to determine the V magnitude associated with 
a given luminosity L, we used the classical magnitude- 
luminosity relation 

V = Mb„i,e - 2.5 log(L/La) - BC + (m - M)v. 

For M13, we adopt a distance modulus (m - M)v = 14.48 , 
following |Cho et al. (2005)| and a bolometric correction BC = 
-0.253, according to the value given by |Girardi et al. (2002")! 



Fig. 12. Same as Fig.^Jfor om rotating reference model Ml. 

for [Fe/H] = -1.5, 7^ = 5000 K and logg = 2. M^.o = 4.76 
according to Cox ( I2000t . 

FigureOemphasises the following points : 

1 . The bump clearly appears for all models in spite of the nu- 
merical noise present at higher magnitudes in models M2 
and M6. When the ju-barrier has been efficiently eroded, 
the bump is less pronounced in the LF. 

2. The predicted LF for models with low degree of mixing 
(namely Ml, our reference model, and M3, with fi^ = 0) 
are similar to the one obtained for the standard model 
MO. In case of model M3, the diffusion coefficient in the 
outer HBS remains always small prior to the bump (i.e. 
Aot ^ 5.10'^ cm-.s"'). This allows preservation of a large /i- 
gradient which is responsible for the clear signature in the 
LF. A similar effect explains the LF for model Ml, where 
transport efficiency is always low in the radiative zone. 

3. In rotating models, the bump occur at higher luminosity 
compared to the standard model MO. 



6.2. Abundances 

In § 16.11 we mentioned the lowering of the mean molecular 
weight barrier left by the 1st DUP in rotating models Ml and 
M6. Figure[21shows the erosion of '^C and '"*N profiles as a re- 
sult of rotational mixing. However, this chemical diff'usion does 
not affect significantly the surface abundances, and we report 
only a minor increase of anti-correlated with a decrease of 
'^C compared to the standard case. 

In Model M6 the higher degree of diff'erential rotation at 
the base of the CE feeds the turbulent shear-induced mixing. 
The outer '"*N plateau is erased and the nitrogen mass fraction 
in the CE is increased relative to the standard case. '^C also 
diff^uses outwards and the peak around 6M = 0.3 is flattened 
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Fig. 13. Same as Fig.[nfor model M6. 

out. Although the mean molecular weight at the depth of deep- 
est penetration of the CE is much lower than in the standard 
model, D,, (which dominates the transport) decreases rapidly 
from 10^ cm^s"' just below the CE down to k IQr' cm^s"' in the 
chemically inhomogeneous regions of the outer HBS (Fig.|S}. 
Mixing is thus confined to a narrow region located just be- 
low the CE, where the chemical profiles are flat. According 
to Charbonnel (1995) and Denissenkov & VandenBerg (2(503), 
diffusion coefficients as large as w 4 • 10*^ cm^s"' are needed to 
connect the HBS with the CE and to modify the surface abun- 
dances. However this configuration is never reached in our self- 
consistent models. 

We reach similar conclusions for the other rotating mod- 
els, since at the bump, the diffusion coefficients have the same 
magnitude, and are too small to affect the surface abundance 
composition. This is clearly illustrated in Fig. 16.21 where we 
present a comparison of the temporal evolution of lithium, car- 
bon, nitrogen and carbon isotopic ratios obtained in our mod- 
els with homogeneous observational data for field stars with 
[Fe/H] € [-2;-!]. While variations associated with the 1st 
DUP are satisfactory and rotating models tend to better agree 
with observations, no further variations are obtained after the 
bump. Let us note that model M6 leads to the destruction of 
lithium already on the MS (because of strong mixing associated 
with the shear in absence of any compensating mechanism) in 
contradiction with the observations. This situation is similar to 
that encountered in Pop 1 stars on the red side of the lithium 
dip. For a complete discussion of this problem, the reader is 
referred to Talon &. Charbonnel f 1998), as well as to Talon & 
Charbonnel ( l2005l and Charbonnel & Talon (l2005i which de- 
scribe how internal gravity waves could help resolve this issue 
in Pop 11 stars. 

Regarding lithium on the RGB the mixing rates associated 
with shear-induced turbulence do not allow the triggering of 



the Li-flash as proposed by Palacios et al. J2001> to consis- 
tently explain the low percentage of lithium rich giants at the 
bump luminosity. Indeed, in this scenario, an initial mixing rate 
of about 10''cm^.s ' is needed in the region of the ^Be peak 
region in order for it to diffuse outwards and decay into 'Li in 
a region where this nuclide will efficiently capture protons and 
give rise to a energetic runaway called the "Li-flash". We never 
get such high mixing rates in the radiative interior of our mod- 
els. 

The present description of the extra-mixing process in RGB 
stars does not allow the models to reproduce the observed abun- 
dance anomalies in upper RGB stars, and does not validate the 
Li-flash scenario for Li-rich RGB. 

7. Comparison with previous works 

In this paper we have presented the first models of rotating 
low-mass stars which take into account rotational transport 
by meridional circulation and shear turbulence coupled 
self-consistently to the structural evolution from the ZAMS 
to the upper RGB. A detailed study of the input physics 
associated with the rotational transport of angular momentum 
and chemicals allowed us to assess the impact of various 
physical ingredients on the extension and magnitude of mixing 
along the giant branch. Let us know compare our predictions 
with others from the literature 

1 . Angular momentum evolution 



With regard to the evolution of the angular velocity pro- 
file, we assumed solid body rotation on the ZAMS in all 
our models and then let angular momentum be transported 
by meridional circulation and shear-induced turbulence. At 
the turn-off, this leads, in all cases, to differential rotation 
in the radiative interior and to a slowly rotating convective 
envelope. Beyond the turn-off, the angular velocity below 
the CE is determined by the assumed rotation law in the 
CE. Both the absolute value o/Qbce <ind the degree of dif- 
ferential rotation (and hence mixing) in the radiative zone 
are maximised after the completion of the 1st DUP in the 
case of a differentially rotating CE ( jcE(i') — est). This sup- 
ports the conclusions anticipated by other authors. In 1979, 
Sweigart & Mengel conjectured that differential rotation in 
the convective envelope of a red giant with homogeneous 
specific angular momentum could be necessary to provide 
enough rotational mixing at the bump. Recently, Chaname 
et al. (12005 > . using a "maximum mixing approach", reach 
the same conclusion from considerations on the global AM 
budget. They also indicate that if Q. has to be described by a 
power law, it is not necessarily with a -2 index (e.g. uniform 
specific angular momentum). As we mentioned in § 13.51 
assuming uniform specific angular momentum is a first ap- 
proximation, and we may expect better estimates of the CE 
rotation regime from direct numerical simulations. 
In their work, Chaname et al. (I2005> nonetheless insist on 
the fact that in order for rotational mixing to reproduce the 
observed abundance anomalies of low-mass Pop 11 RGB 
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Fig. 14. Normalised theoretical luminosity function in V for the models presented in Table 1 and 3. These histograms represent 
the time spent in each magnitude (luminosity) bin as a function of visual magnitude V. The ordinates are normalised to the total 
time spent on the portion of the RGB where V < 16.6. This limit has been chosen in order to have the same normalisation as in 
the observational case presented in the 8''' panel (low right). The observed normalised luminosity function is based on Cho et al. 
(2005). The arrows in panels 3 and 8 indicate the position of the bump. 



stars, they need to assume unrealistic large rotation ve- 
locities at the turn-off. This is a conclusion that we also 
reach in our complete computations. Differential rotation 
in the radiative region separating the HBS from the CE in- 
creases as the CE retreats after the DUP. However, all to- 
gether, the self-consistent evolution of the rotation profile 
for realistic surface velocities at all phases leads to trans- 
port coefficients too low by 3 orders of magnitude com- 
pared to what is expected from parametric studies (Weiss 
et al. \2000l Denissenkov & VandenBere 120051 in order to 
alter the surface chemical composition beyond the bump 
luminosity and reproduce the observed patterns. 
In this work, we did not force the specific angular momen- 
tum to have the same value at all times in the convective en- 
velope , as was done by Denissenkov & Vandenberg (I2003> . 
Under such an assumption the specific angular momentum 
within the CE remains at the same level during the 1st DUP 
as at the turn-off, which results in an increase of the differ- 
ential rotation and mixing below the CE. On the contrary, in 
the models M2 to M6 presented here, when the CE deepens 
during the 1st DUP, it dredges material with lower specific 
angular momentum, and jce (which is the same in each 
mass shell within the CE) drops as can be seen from Fig.Ol 



As a result, the angular velocity and the differential rota- 
tion are lower at the bump. Beyond the completion of the 
1st DUP, jcE continues to evolve slightly due to angular 
momentum transfer with the underlying radiative zone, but 
does not significantly vary anymore. Assuming no varia- 
tion in space of the specific angular momentum in the CE 
after the turn-off appears to be very different from assum- 
ing no variation with time of this same quantity. This latter 
assumption, when combined to strong differential rotation 
in the radiative interior on the MS, is an ad hoc way to pro- 
duce strong differential rotation (and mixing) at the bump. 
Although Denissenkov and collaborators (Denissenkov & 
Tout EUnni (OrOO), Denissenkov & VandenBerg IMHl 
(Dy03)) have also searched for a solution to the RGB 
abundance anomalies problem in terms of rotational mixing 
by meridional circulation and shear turbulence, they reach 
very different conclusions in terms of the evolution of the 
surface abundance pattern on the RGB. As a matter of fact, 
their transport coefficients are very similar to ours in terms 
of shape, but they get at all phases much larger rotation ve- 
locities and differential rotation rates than we do, resulting 
in higher mixing rates. 
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Fig. 15. Evolution of the lithium abundance (e(^^')), the carbon isotopic ratio ''C/'-'C, [C/Fe] and [N/Fe] as a function of the 
luminosity logarithm for the models presented in Table 1. Line encoding is the same as in Fig. 13.51 Black asterisks are actual 
measurements, open upward triangles are lower limits and open downward triangles are upper limits. Observational data from 
Gratton et al. (2000) for field stars in the metallicity range [Fe/H] e [-2; -1]. 



The origin of the differences is difficult to assess although 
several points are certainly critical. First of all the self- 
consistency of the treatment of rotational transport within 
the stellar evolution code, which allows the retro-action of 
AM and chemicals transport on the structure at each evolu- 
tionary step seems to be crucial. Although DTQO solved 
Eq. n using Maeder & Zahn J 19981 1 formalism, this was 
done outside their evolution code in a post-processing way. 
In addition they imposed an angular velocity profile at the 
bump (up to this point their model does not consider any 
transport) so as to reproduce the observed anomalies in the 
globular cluster M92. As a result they obtain very large 
mixing rates able to change the O and Na surface abun- 
dances. As mentioned in §n recent observations in differ- 
ent GCs indicate that the O-Na anti-correlation also exists 
in turn-off stars, with a similar spread as in RGB stars. This 
strongly suggests that this pattern is predominantly of pri- 
mordial origin, and that there is no need for the evolution- 
ary models to reproduce it, at least when considering an av- 
erage GCRGB star''. In our approach, the angular velocity 

^ In the specific case of M13, the extreme O-Na anticorrelation ob- 
served in RGB tip stars could on the other hand be attributed to some 



profile at the bump is not assumed but results from the evo- 
lution (due to structural readjustments and rotational trans- 
port) from the ZAMS. It is by no means a free parameter 
that can be tuned at the bump. 

Another important difference concerns the assumption 
made by DVQ3 on the evolution of the angular velocity at 
the base of the CE (in this paper they do not consider the 
rotational transport of AM, and they apply the rotational 
transport to the chemicals only from the bump on). In 
order to get the "right" Q. profile at the bump, they need to 
impose the constancy of the specific angular momentum 
from the ZAMS up to the RGB tip. We consider that this 
strong assumption is unphysical. 

As a consequence, the work by Denissenkov and collabo- 
rators should not be considered as a proof of the efficiency 
of rotational transport by meridional circulation and 
shear-induced turbulence to modify the surface abundance 
pattern of low-mass RGB stars. It however shows what 
should be the angular momentum distribution at the bump 

extreme mixing whose signature superimposes to the primordial pat- 
tern and exacerbates it. 
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in a rotating star, for shear-induced turbulence to produce 
the expected amount of mixing required in these objects. 

2. Description of the turbulence 



We have investigated the effects of using different prescrip- 
tions for both the horizontal and vertical turbulent diffusion 
coefficients during the RGB evolution. Although some 
differences arise depending on the adopted descriptions 
of turbulence, their effects on the transport of chemical 
species at the bump and beyond are marginal. Concerning 
the choice for v/,, the MPZ04 prescription, predicts a larger 
value which ensures the validity of the shellular rotation 
scheme, a condition that is not always fulfilled when 
using Zahn92 prescription. The MPZ04 prescription also 
(over-)quenches the efficiency of meridional circulation 
and that of the //-currents. Regarding the choice for v,,, 
contrary to the expectations of Denissenkov & Tout 
J2000> . the TZ97 prescription for v,, does not lead to 
efficient mixing throughout the entire radiative zone, nor 
does it contradict the observations when implemented in 
a self-consistent scheme where /i-gradients are taken into 
account. The observations in low-mass RGB stars do not 
allow any discrimination between the MM96 and the TZ97 
prescriptions for v,,, and the use of the former prescription 
in order to prevent any mixing in the outer HBS prior to 
the bump, as advocated by Denissenkov and collaborators, 
is not justified. Let us also recall that the TZ97 prescription 
gives a much better agreement in the case of massive stars. 

3. //-currents and //-gradients 



Several important results obtained from our models 
concern the effect of mean molecular weight gradients 
(//-gradients) on the rotational transport. 

In this paper and in paper I, we have studied the effect of 
the //-gradients on the main sequence. In § 15.31 we have 
shown that if the present low-mass RGB stars were slow 
rotators on the ZAMS, shear-induced turbulence could not 
develop in the radiative interior during the main sequence. 
In this case, and independently of the prescription used for 
the turbulent diffusion coefficients, the //-currents have no 
effects since rotational mixing is negligible. For model M6, 
a fast rotator on the ZAMS undergoing strong braking on 
the main sequence, we reach the same conclusion as in 
paper I when the Zahn92 prescription is used for D/,: //- 
currents play an important role in shaping the turn-off rota- 
tion profile. On the other hand, the use of the MPZ04 pre- 
scription for Dh reverts this conclusion and the //-currents 
are insignificant, as in the case of the slow rotators. 
Beyond the turn-off, the transport erodes the ji-gradients 
in all our rotating models, including those with uniform 
angular velocity in the CE. The //-discontinuity trans- 
lates into a dip in the diffusion coefficient profiles. The 
lesser the mixing, the broader and the more persistent 
this feature. In models with uniform specific angular 
momentum in the CE, this gap is soon filled after the 



1st DUP because turbulent transport is efficient enough 
in this region to smooth the chemical gradients. Despite 
the lowering of the //-barrier, Dtot remains however too 
small to connect the outer HBS with the CE, and the 
surface abundance pattern is not altered prior to the 
bump. Concerning the observational consequences, the 
spreading over of the //-barrier lowers the luminosity 
function height at the bump, but does not erase it. The 
/d-gradients are thus seemingly not entirely responsible 
for the lack of mixing evidence in lower RGB stars, con- 
trary to what was conjectured by SM79 and Charbonnel 
il995\l . A similar conclusion was reached by Chaname et 
al. i2U05l in their "optimised rotational mixing" approach. 

Finally, in all our rotating models, the //-barrier associated 
with the HBS very efficiently prevents any mixing to con- 
nect the Na-rich layers with the outer radiative envelope. 
Thus, contrary to Chaname et al. <2004all2004bll2005l and 
Denissenkov & VandenBerg j2003> the mixing depth does 
not need to be parametrised if the effects of //-currents on 
the mixing are consistently taken into account. 



8. Conclusion 

Our self-consistent approach of rotational mixing associated 
with meridional circulation and shear-induced turbulence leads 
to two major conclusions : (1) this formalism does not pro- 
vide large enough transport coefficients in low-mass, low- 
metallicity RGB stars as required to explain the abundance 
anomalies observed both in the field and in globular clusters; 
(2) it requires differential rotation in the convective envelope^in 
order to obtain non-negligible differential rotation rates, and 
hence mixing rates, in the underlying radiative region. 

These results point toward remaining open questions that 
we would like to bring to light. 

The interplay between convection and rotation in extended 
stellar convective envelopes is still unknown, and the hypoth- 
esis of differential rotation is very attractive when the shear is 
the only process considered to trigger the turbulent transport of 
angular momentum and chemicals. 

In our work, it appears that changing the angular velocity pro- 
file in the convective envelope from uniform to differential, 
increases the degree of mixing in the underlying radiative re- 
gion. This enhancement remains however moderate and does 
not lead to the large diffusion coefficients expected from para- 
metric studies. As the shear-induced turbulence appears not to 
be efficient on its own to reproduce the observed abundance 
anomalies in low-mass giants independently of the rotation law 
in the convective envelope, we are not able at present to make 
any statement concerning this aspect. We nonetheless would 
consider with great interest the use of 3D hydrodynamical di- 
rect simulations to assess the rotation regime within the ex- 
tended convective envelopes of cool giants. 



^ Tliis rotation regime is aciiieved in our case by imposing uniform 
specific angular momentum within this region. 
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In their recent work, Chaname et al. J2005> propose that the 
rotation regime of the convective envelopes may change along 
the evolution, going from solid-body during the MS to differen- 
tial rotation on the RGB. Such a scheme would reconcile rota- 
tion velocities of MS and horizontal branch stars together with 
abundance patterns from the MS to the RGB tip. Here again, 
if another physical process, such as internal gravity waves, is 
able to transport angular momentum also in giants and if its ef- 
ficiency at all evolutionary phases depends on the initial mass, 
the modification of the rotation regime in the convective enve- 
lope during the giant phase could not be necessary anymore. 

The second point concerns the transport mechanisms asso- 
ciated with diff'erential rotation. By now, only the secular shear 
instability has been investigated. This is mainly related to his- 
torical reasons, since Zahn's formalism was at first derived for 
MS stars, where this hydrodynamical instability is dominant. 
The structure of a giant star is however very far from resem- 
bling that of the Sun. We might then expect that specific fea- 
tures such as nuclear burning shell, contracting radiative inte- 
rior and expanding (extended) convective envelope favour the 
triggering of other hydrodynamical instabilities, or physical 
transport mechanisms. Spruit & Knobloch (1984) advocated 
the possibility for the baroclinic instability to be efficient in gi- 
ant stars, but at present, we lack a description in the non-linear 
regime (e.g. in a regime associated with turbulence). This in- 
stability also depends on the degree of differential rotation but 
is not as sensitive to //-gradients as the secular shear, and could 
thus complement the effect of secular shear to increase the de- 
gree of mixing in regions with large yu-gradients. Other insta- 
bilities such as the Goldreich-Schubert-Fricke (GSF) and the 
Solberg-H0iland instabilities could also become non-negligible 
during the RGB phase. 

Last but not least, and in connection with the presence 
of the so-called "super Li-rich giants" at the RGB bump 
(Charbonnel & Balachandran 120001 and references therein), 
Palacios et al. J200H proposed that the structural and nuclear 
response of the star to the rotation-induced mixing could cause 
an increase of ZJtot as required to explain the abundance anoma- 
lies. This scenario involves an important release of energy due 
to ^Li burning in the external wing of the HBS (the so-called 
"Li flash"). This opens a new field of investigation concerning 
the reaction of meridional circulation and of the various insta- 
bility to a major and local release of nuclear energy. 

Our results clearly show the lack of success of secular 
shear alone to trigger large enough amount of extra-mixing, 
but does not rule out rotation-induced mixing being responsi- 
ble for abundance anomalies in low-mass giants. As a matter of 
fact, rotation remains at present the best candidate for the trig- 
gering of the extra-mixing process in these objects. Our study 
points toward the need for more physics, implying a revision 
of the coupling between rotation and convection, and the treat- 
ment of other physical processes likely to contribute to trans- 
port angular momentum, including hydrodynamical instabili- 
ties neglected before. 
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Appendix A: Detailed expressions for the 
components £n and in Eq.|51 

Maeder & Zahn (I1998I I derived a generalised expression of the 
meridional circulation velocity [/,., taking into account the ef- 
fects of yu-gradients, non-stationarity and of a general equation 
of state, that was revised and corrected in Denissenkov & Tout 



(I200(H . 

Here, we reproduce the detailed expressions of the and 
Efj terms appearing in Eq.|2] Eauations lA. 1 l and lA. 3l correspond 
to corrected versions of Eq. 4.30 and Eq. 4.42 found in Maeder 
& Zahn a998j . 
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In these expressions we have used the following notations: 



dlnT 

is the temperature scale height, 

X 4acT^ 



K — 

pCp 3p-KCp 

is the thermal diffusivity and x is the thermal conductivity, e 
and E^'"^ are the means over an isobar of the nuclear and grav- 
itational energy respectively, and 

represents the nuclear fraction of energy. e„, -L{r) I M{r) and 
p,n - M{f)l{Anr' IJ) are the mean total energy and the mean 
density inside a sphere of radius r. 



^/51n^\ 



' d\nx \ 

are logarithmic derivatives of the thermal conductivity x , and 



dlnE\ ( 51ne 



are the logarithmic derivatives of the nuclear energy. 

d\np \ ^_ ld\np\ 

are derived from the equation of state using the notation by 
Kippenhahn & Weigert ( I1994> . 

The term Ea. lA.ll is of order 2 in Q- whereas the 

other terms are of order 1, and is thus neglected. 



